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We approach the biexciton Schrodinger equation not through the free-carrier basis as usually 
done, but through the free-exciton basis, exciton-exciton interactions being treated according to the 
recently developed composite boson many-body formalism which allows an exact handling of carrier 
exchange between excitons, as induced by the Pauli exclusion principle. We numerically solve the 
resulting biexciton Schrodinger equation with the exciton levels restricted to the ground state and 
we derive the biexciton ground state as well as the bound and unbound excited states as a function 
of hole-to-electron mass ratio. The biexciton ground-state energy we find, agrees reasonably well 
with variational results. Next, we use the obtained biexciton wave functions to calculate optical 
absorption in the presence of a dilute exciton gas in quantum well. We find a small asymmetric 
peak with a characteristic low-energy tail, identified with the biexciton ground state, and a set of 
large peaks associated with biexciton unbound states, i.e., exciton-exciton scattering states. 
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I. INTRODUCTION 

The existence of composite particles in semiconduc- 
tors has been predicted long agc-k Bound states made 
of conduction band electrons and valence band holes re- 
sult from the Coulomb attraction between these carriers. 
To name the simplest ones, these bound states are exci- 
tons (X) made of one electron plus one hole, trions (X T ) 
made of two electrons plus one hole, or two holes plus one 
electron, and biexcitons (XX) made of two electrons plus 
two holes. More exotic composite objects made of a large 
number of correlated fcrmion pairs, called "electron-hole 
droplets" , also exist, with a carrier density far larger than 
the one in which excitons are formed. 

The simplest composite particle, the exciton, made 
of one conduction electron and one valence hole, is 
very similar to an Hydrogen atom, if we neglect in- 
terband Coulomb processes 2 -. Exciton has bound and 
unbound (scattering) states which can be analytically 
determined 3 -. Exciton bound states appear as large nar- 
row peaks in the photon absorption spectrum, the peak 
intensity depending on the so-called "exciton oscillator 
strength". The reason bound exciton peaks are easy to 
observe is twofold: first, when a plane- wave photon with 
momentum Q transforms into a bound exciton, the cou- 
pling is quite good because the center-of-mass motion of 
the bound exciton also is a plane wave with same mo- 
mentum Q. Second, excitons, made of an even number 
of fermions, have a bosonic nature; so, they can be piled 
up all at the same energy, the absorption peak intensity 
increasing linearly with the density of excitons already 
present in the sample. 

Observation of composite particles like trions is more 
complex. It has been hampered for quite a long time, 
partly due to the small binding energies that trions have 
in bulk samples, at best one order of magnitude smaller 
than the exciton binding energy. Such binding energies 



are smaller than usual exciton line width; so, trion peaks 
fall on the side of exciton lines. This small trion binding 
energy can be physically understood by seeing the trion 
as an electron or a hole bound to an exciton. The effec- 
tive attraction then is dipole-like which makes interaction 
between exciton and free carrier much weaker than be- 
tween elementary charges. A clear signature of trions has 
been obtained recently only in semiconductor quantum 
wells^— , the reduction of dimensionality increasing all 
binding energies as seen from the exciton energy which 



goes from i?^ D ' in 3D to R 



? (2D) 
X 



4i?^ D) in 2D, and to 



infinity in ID. 

What makes bound trions hard to observe has also to 
be traced back to their oscillator strength which is one 
trion volume divided by one sample volume smaller than 
the exciton oscillator strength^. This drastic reduction 
factor can be physically understood as the probability for 
a photocreatcd exciton to localize over a trion volume, a 
free carrier initially spread over the whole sample. As 
a result, the trion peak commonly observed in heavily 
doped samples in which a large electron density exists, 
should not be interpreted as a signature of elementary 
trion, but rather as an exciton interacting in a coherent 
way with all the electrons present in the sample. Such 
a many-body effect is singular and leads to a broad ab- 
sorption line, as experimentally showrA 

Mathematically, the derivation of trion eigenstates 
amounts to solving a three-body problem which has no 
known analytical solution. Attempts to tackle such a 
problem inevitably rely on some truncation scheme in 
addition to heavy numerics in order to possibly obtain 
satisfactory results. Recently, we showed how, using a 
physically relevant viewpoint, we can approach one trion 
as an exciton interacting with an electron-^. We have 
constructed a Schrodinger equation for trion using the 
electron-exciton basis and solved it by restricting this 
basis to the first few low- lying exciton levels, — which is 
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reasonable since the energy scale for excitation of the ex- 
citon internal motion, of the order of one Rydberg, is 
much larger than the other energy scales. The resulting 
trion binding energies we find in 2D and 3D agree reason- 
ably well with the most accurate variational results. One 
important advantage of this approach is to allow reach- 
ing the trion ground and excited states on equal footing, 
these excited states being out of reach from standard 
variational methods. 

Following Lampert's prediction^, an even more com- 
plex composite particle, the biexciton, has been observed 
in bulk materials such as CuCl [Ref.[n[, Cu 2 [Ref.[lJ], 
and AgBr [Ref . GUH . More recently, the biexciton bind- 
ing energy has been measured in GaAs quantum welU 5 - 
and found to be one order of magnitude larger than in 
bulk samples, a result supported by calculations done 
one year later—. Since then, other aspects of biexcitons 
in confined structures, such as optical enhancement in 
biexciton formatio n 17 ' 18 and the influence of dimension- 
ality on the biexciton bindin g 15 ' 19 ' 20 , have been studied. 
Biexcitons in quantum wires have also been reported^. 

In view of the successful application of the compos- 
ite boson many-body formalism to trioni -, we, in this 
work, go on along the same line to tackle biexciton. The 
biexciton problem a priori is an even more complex four- 
body problem, with two electrons (ei,ea) and two holes 
(hi, I12) involved. The idea is to start with two electron- 
hole pairs bound into two excitons by the strong electron- 
hole Coulomb attraction. The exciton-exciton attraction, 
although quite weak since it essentially is dipole-like, al- 
lows two free excitons to form a molecule with a binding 
energy substantially smaller than the exciton binding en- 
ergy. To approach a system made of two electrons plus 
two holes, the exciton basis is physically quite appeal- 
ing because the strong exciton binding energy is then 
included into the problem at the zeroth order. We are 
left with solving a Schrodinger equation for the weaker 
biexciton binding energy. In this approach, the four-body 
system is pictured as two interacting excitons: one exci- 
ton is made of the (ei, h\) pair, while the other is made of 
(e2, /12) pair, these two excitons however exchanging their 
carriers to be possibly made of (ei, /12) and [e-2, hi). Such 
a two-exciton picture could be thought, at first sight, to 
lead to an easy problem because of the weak exciton- 
exciton attraction compared with the strong electron- 
hole attraction. However, this weak attraction is the one 
responsible for two excitons to be bound into a molecule. 
So, in order to reach bound states and find the associated 
poles, this exciton-exciton interaction has to be treated 
in an exact way. 

With this goal in mind, we here construct a biexciton 
Schrodinger equation in terms of the exciton basis us- 
ing the recently developed composite boson many-body 
theory^. In much the same spirit as Feynman diagrams 
for elementary particles, this theory takes advantage of 
"shiva diagrams" to visually identify many-body effects 
involved among composite particles. It moreover enables 
treating exactly carrier exchange which results from the 



indistinguisability of the fermionic components of these 
composite particles. By restricting the exciton levels to 
the ground state only, it becomes possible to numeri- 
cally solve the biexciton Schrodinger equation quite eas- 
ily. The values we obtain for the biexciton ground state 
energies in 2D and 3D are in good agreement with varia- 
tional results. One important advantage of the procedure 
is that the biexciton Schrodinger equation can be cast 
into a generalized eigenvalue problem; so, we can reach 
bound and unbound excited states at once, with a single 
matrix diagonalization. 

In a second step, we use the obtained biexciton relative 
motion wave functions to calculate the photon absorp- 
tion spectrum in quantum wells, assumed to be exact 2D 
systems. Instead of considering biexciton as generated 
through two-photon absorption 2 ^—, we here study one 
photocreatcd exciton interacting with a dilute exciton 
gas. Using similar arguments as those we used for bound 
trion, we find that the biexciton oscillator strength is one 
biexciton volume divided by one sample volume smaller 
than the exciton oscillator strength. This would make 
observing the biexciton line very difficult. However, biex- 
citons, like excitons, are boson-like particles: They can 
thus be packed up all at the same energy level. As a re- 
sult, the biexciton absorption line increases linearly with 
exciton density provided that the density is low enough to 
possibly neglect many-body effects between the photocre- 
ated exciton and the free excitons present in the sam- 
ple. The calculated photon absorption spectrum shows a 
small peak, with a characteristic low-energy tail, originat- 
ing from the biexciton molecular state, and large peaks 
centered on the exciton ground levels, which are asso- 
ciated with exciton-exciton scattering states. Both, the 
bound and unbound biexciton peak intensities decrease 
when the temperature increases. It can also be shown 
that the intensity of the absorption line for one biexci- 
ton made from a photocreatcd exciton and an exciton of 
the exciton gas, increases linearly with photon number 
and exciton density. This is in contrast to the biexciton 
absorption line associated with two-photon absorption 
which increases quadratically with photon number and 
thus becomes dominant at high laser intensity. 

The present paper is organized as follows: 

In Sec. HH we briefly discuss the relation which exists 
between biexciton written in the free-carrier basis, and 
biexciton written in the exciton basis. We also intro- 
duce the four commutators necessary to properly handle 
many-body effects involving composite excitons. 

In Sec. IIII1 we study triplet biexciton states made of 
same-spin electrons and same-spin holes and we derive 
the corresponding Schrodinger equation. 

In Sec. HVl we study singlet and triplet biexciton states 
made of opposite-spin electrons and opposite-spin holes. 
These biexciton states are first constructed in terms of 
two free electrons plus two free holes, and then in terms 
of two free excitons, in order to reveal important par- 
ity relations. We then concentrate on singlet biexciton 
states with center-of-mass momentum equal to zero and 
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FIG. 1: (a) \h( n 3 . ) = A(™ j . ) for hole exchange, the exci- 
tons m and i having the same electron, (b) Pauli scattering 

^ e ( m \ ) = ^( r i ) ^ or e l ec t ron exchange, the excitons m and 
i having the same hole. 



we restrict the exciton levels to the ground state. This 
nicely reduces the biexciton Schrddinger equation to a 
ID integral equation. 

In Sec. [Vj we numerically solve this ID integral equa- 
tion to obtain the biexciton binding energies for the 
ground and excited states as a function of hole-to- 
electron mass ratio. We also show the biexciton relative 
motion wave functions for the bound state as well as for a 
few unbound states. Finally, we use these wave functions 
to calculate the photon absorption spectrum in the pres- 
ence of a dilute exciton gas for various low temperatures. 

In the last section, we conclude. 



II. BIEXCITON ON THE EXCITON BASIS 

We consider a system made of two electrons and two 
holes in a semiconductor: the electrons carry a spin 
s = ±1/2 while the holes carry an angular momentum m 
that we will also call spin. In bulk semiconductors, the 
hole angular momentum can be m = (±3/2, ±1/2), while 
in narrow quantum wells, it reduces to to = ±3/2 due to 
the heavy-light hole energy splitting induced by the well 
confinement. For simplicity, here we shall neglect the 
role of light holes and consider heavy holes only. Fur- 
thermore, we neglect the warping of the semiconductor 
valence band and approximate it by a spherical, parabolic 



FIG. 2: Pauli scattering \( ,' , 1 for carrier ex- 
change between an exciton i = (fj,Q) and an exciton j = 
(i>j, — Q) (see Eq. (|B.4|) '). The exciton (i/j, Q) is a linear com- 
bination of electron-hole pair (k + a e Q, k + cenQ) where 
a e = 1 — ecu = m e / (m e + m h ) (see Eq. (|A.7|) 1. 



band. The usual basis for such a biexciton system is then 
made of states with two free electrons and two free holes 
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To transform this free-carrier basis into an exciton basis, 
we make use of the relations which exist between free 
electron-hole pair creation operators and exciton creation 
operators, namely, 



l]Si7Tli 



bl_ „.(k k k R |i), (2) 
i|k e k h ), (3) 



J2 a L,s 
k c k h 



where \i) denotes the i exciton state. Using Eq. ([3]), 
we can rewrite the two-free-electron-hole pair states of 
Eq. |T]) in terms of exciton states as 



(4) 



with (si,Sj) — (s±, S2) and (mi,mj) = (m^ma). 
Note that the basis made of [(si, mi); (S2, ^2)] and 
[(si, TO2); (S2, rrii)] are equally valid. This means that, for 
si =/= S2 and mi ^ m2, the basis can be made either of two 
bright excitons, (-1/2,3/2) and (1/2,-3/2), or of two 
dark excitons, (1/2,3/2) and (-1/2,-3/2), the bright 
exciton basis however being more convenient for prob- 
lems dealing with photons. Note that bright and dark 
excitons are degenerate if we neglect intcrband Coulomb 
processes. 

While the great advantage of the exciton basis is to 
contain part of the electron-hole interaction, actually the 
strong part leading to exciton bound states, its main dis- 
advantage is to be overcomplete; as a direct consequence, 
this basis is not orthogonal. It is possible to overcome the 
difficulties induced by the non-orthogonality of the ex- 
citon basis through the commutation technique recently 
developed for composite boson many-body effects^. The 
keys of this formalism rely on just four commutators be- 
tween exciton operators B\: Fermion exchanges follow 
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from two commutators which read, in the absence of spin 
degrees of freedom, as 
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[M^') + M^M- (6) 



D m i is called "deviation-from-boson" operator because, 
without it, B\ would reduce to an elementary boson op- 
erator. The Pauli scattering Xh ( n J - ) corresponds to a 
hole exchange between cxcitons the excitons m and 

i having the same electron, as defined in Eq. (|B.1[) of the 
appendix and shown in the diagram of Fig. 1(a) In the 
same way, the Pauli scattering A e ( n °. J corresponds to 
an electron exchange, the excitons m and i having the 
same hole, as defined in Eq. (|B.2[) and shown in the dia- 
" Xh ( "j { ) , which is equal to X e (™{), 
ritten as A ( ^ \ ) for simplicity. 
The other two commutators that handle fermion- 
fermion interactions, are 



gram of Fig. 1(b) 
is often wr 



H,B}] = EiB} + Vf, (7) 
tf.-B) = E^ ir (™i)^- (8) 



The "creation potential" V} generates, through (JSJ), 
the direct Coulomb scattering £ dlr ( n 3 ) of excitons i 
and j. It consists of four Coulomb processes between 
the fermionic components of these two excitons: one 
electron-electron repulsion, one hole-hole repulsion, and 
two electron-hole attractions, as shown in the diagram of 
Fig. [3] The precise expression of the hole-hole part of 
this scattering is given in Eq. (|B.5[) of the appendix. 

These four commutators are used to get the biexciton 
Schrodinger equations derived in the next sections. To 
include the electron and hole degrees of freedom, let us 
focus on the two relevant cases: 

(a) Carriers with same spins, s% = S2 and mi = rri2- 
This corresponds to triplet states for the electron part 
(S e = ~i-,S' = ±1) and triplet-like states for the hole 
part (Sh = 3,5^ = ±3). The associated orbital wave 
functions then have to be odd with respect to exchange 
of the two electrons or the two holes in order to fulfill the 
Pauli exclusion principle. 

(b) Carriers with opposite spins, s± = — S2 and mi = 
—mi- The resulting spin configuration then depends 
on the way electrons and holes are linearly combined: 
We can either have a triplet state for the electron part 
(S e = 1, SI = 0) and a triplet-like state for the hole part 
(Sh = 3, S^ = 0), or a singlet state for the electron part 
(S e = 0, Sg = 0) and a singlet-like state for the hole part 
(Sh = 0, S% = 0). This case thus requires a more careful 
derivation since the orbital wave function for triplet state 
must be odd as in the case (a), but even for the singlet 
configuration. The biexciton ground state belongs to the 
set of singlet states. 
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FIG. 3: Direct Coulomb scattering £ "(™ J .) between exci- 
ton i and exciton j. The "out" exciton m is made with the 
same electron-hole pair as the i exciton. Similarly for exci- 
tons n and j. This exciton-exciton scattering consists of four 
terms: two repulsive interactions between electrons and be- 
tween holes, and two attractive interactions between electron 
and hole. 



III. BIEXCITON MADE OF ELECTRON-HOLE 
PAIRS WITH SAME SPIN si = s 2 ,mi = m 2 

Let us start with triplet biexcitons made of same-spin 
electrons and same-spin holes and drop the spin indices 
to make notations of this section lighter. We look for the 
biexciton eigenstatcs 

(ir = (9) 

in the two-free-exciton basis \ij) = BjB!:\v), namely, 



Y,$ ] B\B]\v). (10) 



Since B\B'j = B^Bl, we can replace the above prefactor 
by (4>if + ^ji)/ 2 - Thc biexciton state then appears as 
in Eq. (fT0|) but with the symmetry condition <fy^V = (/>£■ . 

Equations ([7]) and © allow us to rewrite the biexciton 
Schrodinger equation © as 
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with E rs = E r + E s . In the standard case, i.e., when the 
basis is made of orthogonal states, thc above equation 
forces the bracket to be zero. The situation is more subtle 
with the exciton basis because, due to Eqs. ([5|) and (|6|), 
the scalar product of two-exciton states reads as 
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FIG. 4: Part of the direct Coulomb scattering 
£ dir ( ("' q')' - q) ) between a i = (z^Q) exciton and 
a J = Q) exciton coming from hole-hole repulsion 

(second diagram of Fig. O • 



By projecting Eq. (fTTl) onto 



we then find, since 
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= (E mn - 4 

) is defined as 
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where f M ( ™ j 



e ( " ) (™i) = ^ dir (™i)-r(^i)-A(^i)(^-£, ) ). 

(14) 

with = E, s A(^^)e dh '(^") being "in" 

exchange- Coulomb scattering with Coulomb interactions 
taking place between the "in" exciton pair i.e., be- 

fore hole exchange (see the diagram of Fig. [5J. 

£ (mi) can a PP ear as an effective Coulomb scatter- 
ing between two excitons although it also depends on the 
biexciton energy £ n . The first two contributions to this 
effective scattering correspond to the standard combina- 
tion of Coulomb processes appearing, for example, in the 
time evolution of two excitons. The third term is more 
interesting because it only comes from the Pauli exclusion 
principle. Since the associated Pauli scattering A ( ) 
is dimensionless, it must go along with an energy to pro- 
duce an energy-like scattering, this energy actually being 
an energy difference in order to be gap independent. 

We now introduce a similar "out" exchange- 
Coulomb scattering £ out ( " j ) , defined as £ out ( nj ) = 
J2rs £ du (mr)A( s "0 ; Coulomb interactions taking place 
between the "out" exciton pair (m, n) once exchange 
has occurred. It is easy to check that £ m ( " \ ' ) — 

[£° Ut (im)]*- Moreover, the "in" and "out" exchange- 
Coulomb scatterings are not independent: their differ- 
ence is related to Pauli scattering via 

r ( i. i ) - c out (zi) = ( E mn - By) a ( » i ) , (is) 

as easy to recover from calculating the matrix element 
{v\B m B n H B\B^\v) with H acting either on the right or 
on the left. 

By inserting Eq. (|15[) into Eq. (|14[) . we can symmetrize 
the biexciton Schrodinger equation (fT3"|) for triplet states 
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FIG. 5: "In" exchange-Coulomb scattering f in ( ^ 3 ) between 
exciton i and exciton j. The exciton pair exchanges their 
holes after Coulomb interactions, while the excitons m and i 
keep the same electron. 
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FIG. 6: Part of the "in" exchange-Coulomb scattering 
£ m { (v' q') (t 'q*^') between a i = (^i,Q) exciton and a 
j — (fj, — Q) exciton, coming from hole-hole repulsion (first 
diagram of figure [SJ. The exciton pair exchanges their holes 
after Coulomb interactions (see Eq. ()B.15|l ). the excitons m 
and i keeping the same electron. 



a.s 



= (E mn £ v )^± + ( ™ { ) </> 

ij 



(')) 
ij i 



(16) 



the effective Coulomb scattering, defined as 
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now being symmetrical with respect to the "in" and 
"out" exciton pairs and (m,n). 
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IV. BIEXCITON MADE OF ELECTRON-HOLE 
PAIRS WITH OPPOSITE SPINS 

Si = — 82, m\ = — m,2 



The biexciton ground state belongs to the set of singlet 
states (S e = Sh = 0). 



When the electron spins and hole spins are opposite, 
the derivation of biexciton singlet and triplet states with 
$e = &h = ^ requires a more careful analysis of the 
parity condition induced by the free-carrier fermionic na- 
ture. To this end, we first construct biexciton eigenstates 
in the free-carrier basis and derive the parity condition 
imposed by the Pauli exclusion principle for singlet and 
triplet wave functions. We then use Eq. to rewrite 
the biexciton eigenstates in terms of exciton operators 
and rederive the parity condition in the exciton basis. 
Finally, we write down the biexciton Schrodinger equa- 
tion. 



A. Free-carrier basis 

The biexciton eigenstates for opposite electron spins 
and opposite hole spins can be written on free-carrier 
states as, 
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This writing covers singlet state, (S e = Sh = 0), as well 
as "triplet" state, (S e = I, Sh = 3), for which we have a 
sum instead of a difference of pair states — the biexciton 
triplet state S z = being degenerate with respect to the 
ones constructed on (si = 82, m\ = 1112) in Sec. IIIII 
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where ^J/^^ = Vl.k^.k; 
follows the parity condition 
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If we do the same for the hole part, we end up with 
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has the expected parity condition, namely 

= (- 1 ) Sh ^k e ,k^,k^,k ft - 
(22) 



k (n,s e ,s h ) _ (_ u s e An,Se,s h ) _ ( ^s kA (v,Se,s h ) 
- \ l ) <Pk' k e ,k h ,ki 



'k e ,k' kfc.k' 



B. Exciton basis 

To write the biexciton on the exciton basis, we use 
Eq. ([3]) to rewrite the electron-hole state of Eq. (|2"Tj) in 
terms of bright exciton operators. We find 



i4:k)=E 



Av,S c ,S h ) D f of 



Bl^B^v), (23) 



where the prefactor 



,(n,s e ,s h ) 

'ij ' 



defined as 



$ AA) = E E <<l W0Xk A >^£ k , , 

k e ,k' h k^k,, 

(24) 

is the biexciton wave function "in the exciton basis" . 

Using (kfckelPeP/j) = 6 kcP J khPh and the exciton clo- 
sure relation, |*)(*| = it is easy to show that from 
Eqs. l[2"4"j). (jB.lj) . and (|B.2|) . the parity conditions for elec- 
tron exchange or hole exchange read as 

*g£f- A) = (-i) s -£M^>S A,5h) 

= (-l) Sfc E^( n "i)^' Se ' Sh) - (25) 

By noting that ^h{mi) = ^ e ( n i ) » these two equations 
give the parity condition for exciton exchange as 

<$> s » s ») = ( 26) 

We now turn to the biexciton Schrodinger equation in 
the exciton basis. We find, from the commutators ([7]) 
and ©, 

= (H-4 S '' Sh) )l4!k> (27) 

= £{(^.-4 fl - ,sfc) )^ ,s " A) 

+ E^ dir (^')^' Se ' Sh) }^,-i<iK'>. 

jj 

As (v\B n ,iB m ,^iBl _ 1 Bl 1 \v) reduces to 5 m>r &n,a for ex- 
citons made of carriers with different spins, the projection 
of this equation onto {v\B n ^iB m ^i simply gives 

= (E mn ^. A) + ^ ^ (: i ) 

y 

(28) 

Note that, if we instead project it onto {v\B m ^B n} -i, we 
end up with the same equation but with m and n inter- 
changed. So, the resulting Schrodinger equation (|28p has 
to be solved self-consistcntly with the parity conditions 
(|25p . — which is not convenient numerically 
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Moreover, the Schrodinger equation (|28|) for spin 
triplet states (S e = 1, Sh = 3) does not readily reduce to 
Eq. (p~3|) . whereas the Schrodinger equation must be the 
same for all triplet states since they are degenerate. To 
possibly relate these two Schrodinger equations and also 
avoid handling the parity conditions (|25[) . we can intro- 



duce two new sets of functions L p^.' ,hS, " Sh ^ with c = (e, h) 



and rewrite 



MS 



ill It 



(29) 

so that the parity condition ([25)1 is automatically ful- 
filled whatever ( p£r' , ?'' Sc '' s, ' i )_ gy inserting the above equa- 
tion into Eq. ()28l) . we find a Schrodinger equation for 



(pte"' 8 " 8 "). It reads 







+ E[^ dir (:^)+(-l) S ^° Ut (:D (30) 



+(-i) 5 -(^-4 s " flfc) ) Ae («') 



tout n j 



If we now use Eq. (fl5|) to rewrite £' 
f 11 !" 3 ), the above equation also reads 



(c; v ,S e ,S h ) 



) in terms of 



n - (E -£(Se,S h )) (cyr,,S e! S h ) , t(c;77) fnj) 
u ~~ l^mn °?7 J'rmn ' / 4 Ssym \rni) 



9, 



(31) 

where the effective scattering now has a symmetrical form 
with respect to "in" and "out" states 



<:(c;f) ( n j- 

Ssym V m i . 



= € 



dir f n j 

m i , 



(-iy 



■out f n j 

m i . 



+A C 



(32) 



and zero otherwise. 



It is then easy to see that, for spin triplet state, Sh = 3 
and c = h, the Schrodinger equations (|3"Tj) and (TT3")) arc 
indeed identical. 

Equations (fT5|) and (|2"g)l . in the absence of Coulomb 
scatterings £ dlr ( ™ •? ! ) and Pauli scatterings A 

would lead to ^ = 1 for E mn 
The biexciton would then reduce to two free excitons as 
expected. Interactions between excitons, through both 
Coulomb and Pauli scatterings, produce a difference be- 
tween the biexciton energy and the energy of two free 
excitons. 

By comparing the Schrodinger equations (fT6|) and (fl7|) 
for same-spin carriers with the Schrodinger equations 
(f3"Tj) and (|52")l for opposite-spin carriers, we see that there 
is a sign change in front of the exchange part of the ef- 
fective scattering (|32|) . depending on singlet and triplet 
states. This sign change which results from the Pauli 
exclusion principle, is the reason why two excitons in a 
singlet state (S e = Sh = 0) can bind together into a 



molecular state. This can be understood by considering 
energies close to the energy 2e Vo of two ground state exci- 
tons labeled by (vq, Q = 0), which is the expected biexci- 
ton energy for temperature much smaller than the energy 
for exciton internal motion excitation. So, the biexciton 
binding energy 8 n = 2e„ — £^ h ' is expected to be far 
smaller than the difference between exciton binding ener- 
gies e Vi —s Uo . We then note that the direct Coulomb scat- 
tering £ dlr (o o) i s e q ua l to zero; so, £ dlr essentially van- 
ishes for small momentum transfer (see Eq. (|B.10p ). By 
contrast, the "in" exchange-Coulomb scattering is given 
by2L2£ 



; m (88) = 



- 8?r- 



315tt 3 \ (ax\ 2 u 



512 

267r fax^ 3 
3~ \~T 



DO 1 

V L J Rx 



R 



(3D) 
X 



in 2D 



in 3D 



out / 



with f (§§) equal to £ 



(33) 

according to Eq. (|15l) . 

As a result, the sum of the "in" and "out" exchange- 
Coulomb scatterings, £ in (oo) renders the 
effective scattering ^sym overall largely negative for small 
exciton momenta, while the third term of Eq. (|32[) , of the 
order of the biexciton binding energy, is small. This large 
negative effective scattering allows bound-state solutions 
to the Schrodinger equation (|5Tj) . It is then crucial to 
treat the exchange-Coulomb interactions adequately if 
one aims at getting reliable results for the ground and 
excited states. Detailed discussions about the depen- 
dence of Coulomb scatterings on electron-to-hole mass 
ratio and on relative motion momentum of the exciton 
pair can be found in Refs. [29| and l30l. For completeness, 
in | Appendix II[ we have rederived the various scatterings 
appearing in the biexciton Schrodinger equations. 



Equation 



for i — (vi , k + K/2) and j = , — k - 



K/2), K being the biexciton center-of-mass momentum 
and k the relative motion momentum between the ex- 
citon pair, allows us to rewrite the biexciton operator 
in terms of two bright excitons as (see also Eq. (29) of 
Ref.ED) 



E 



Av,s*,s h ) 



[Vii u j)^t u \ SL +K/2;-l^l j ,-\si+-K/2;V 
Vi Vj ;k 

(34) 

Since all A and £ scatterings do not depend on the center- 
of-mass momentum K, we can, without any loss of gen- 
erality, set K = from now on. Such a biexciton is then 
made of two excitons with opposite momenta. 

Since the parity condition (|25[) is difficult to numer- 
ically implement in the Schrodinger equation (|28[) for 
^n,s c ,s h ) we instead solve a somewhat more com- 

plicated equation (pH]) in which the parity condition (|2"9")l 
is already enforced, the function ^^-.S^Sh.) ^ ^ ^ 
Eq. (|2T)|) being a free parameter. 

Let us from now on focus on the biexciton singlet 
state (S e = Sh — 0) and consider c = h without any 
loss of generality By restricting the exciton level to the 
ground state vq and by setting ( ^' l ' I? ' ' ) ( l / 0) j/ ) = tp^\ 



it 
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the Schrodingcr equation (|3Tj) for the singlet state then 
reduces to 



x ( (Ko,-k) (fo-k'A .An) 
* \ (v M) K,k') ) *V 



M 



X 



:*4 ,) +Ef( ( SS )( S^ ) )^ (35) 



where Mx — TO e + W/j and 



(i/ ,— k) (i/ ,— k') 
(i/ ,k) (^o,k') 

1 r 

f 2 



e dir ( 



(i/ ,-k) 0o,-k') 
(i/ ,k) (i/ ,k') 



(i/ ,-k) 0o,-k') 
(y ,k) (i/ ,k') 



r ut ( 



Oo--k) {v ,— k 
(i/ ,k) (fo,k 



'?) 



-A 



fo.-k) (* ,-k')\ ( k ' + k2 ) 



(i>o,k) («Oik') 



71/ 



(36) 



The scatterings A and £ depend on (k, k') through 
|k|, |k'| and the angle 6>kk' between k and k', the explicit 
values of these scatterings being given in |Appendix II 



Since the biexciton singlet state has a zero angular mo- 
mentum, the (p^ function depends on |k| = k only. So, 
we can first average the various scatterings over the #k, k' • 
Let us call A(fc, k') and £(k, k') these averaged quantities. 
We then end up with a ID integral equation for the func- 
tion ipy' . It reads as 



M x 



(37) 



This equation is numerically solved to get the binding 
energies of the biexciton ground and excited states for 
various hole-to-electron mass ratios. We can also get the 

ifip functions, which are related to the biexciton wave 
functions with proper symmetry via Eq. (|29p . The 
excited states are expected to mainly come from vibra- 
tional modes. To reach rotational modes, it is necessary 
to include p and d exciton levels. As a consequence, the 
A and £ scatterings, as well as the biexciton wave func- 
tions <j)jj , would get an angular dependence. The precise 
treatment of the angular dependence in these scatterings 
is rather complex and definitely beyond the scope of the 
present work. A relatively simple, yet nontrivial, biexci- 
ton state follows from just considering two ground-state 
excitons with a nonzero relative-motion angular momen- 
tum. 



V. RESULTS AND DISCUSSIONS 

To solve the Schrodinger equation (|3"T1) for the biex- 
citon binding energy S^, we use the normalized exciton 
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FIG. 7: (Color online) Binding energies of the 2D biexciton 
ground state(GS) and the first three bound excited states in 

(3D) 



p(2D) 

mass ratio mj,/m e 
minimum for m e 



= AR^> unit, as a function of the hole-to-electron 
The ground state binding energy has a 
m h . 
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FIG. 8: (Color online) Same as Fig.[7]for 3D; the energy unit 
now is /4 3D) . The curves are qualitatively similar to the 2D 
curves, though the binding energies are significantly smaller. 



ground-state wave functions (p^o) in 2D and 3D 



<p|ls) (2D) 
<p|ls} (3D) 



fax' 

V L J (l + a|p 2 /4) 3 / 2 ' 

/0£\ 3 / 2 80T 

\l) (l + a^p 2 ) 2 ' 



(38) 
(39) 



L is the sample size, ax = fi^sc/^xe 2 is the 3D ex- 
citon Bohr radius, with /i^- 1 = m~ x + , and e sc is 
the static semiconductor dielectric constant, one order of 
magnitude larger in semiconductor samples than in vac- 
uum. 

The Schrodingcr equation (|3~T|) can be cast into a gen- 
eralized eigenvalue problem, with matrices spanned by 
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the k momentum. To solve it, we sample the k = |k| 
values with 150 mesh points in 2D and 100 in 3D, ac- 
cording to hi = uf where the Mi's are equally distributed, 
thereby allowing for more sampling in the small k region. 
The upper cutoff /c max (in a^ 1 unit) is taken to be 10 in 
3D but 20 in 2D because the exciton wave function has 
a larger radial extension (in k space) in 2D than in 3D 
and also because the Coulomb interaction Vq decreases 
as 1/q in 2D, i.e., more slowly than the 1/q 2 dependence 
it has in 3D (see Eq. (|X6lO . 



A. Biexciton binding energies for ground and 
excited states 

Solid curves in Figs. ([7]) and ([8]) show the biexciton 
binding energies in 2D and 3D as a function of hole-to- 
electron mass ratio rrih/m e . The results are expressed in 
terms of the corresponding effective Rydbcrgs, namely 
i4 2D) and 4 3D) , with i4 2D) = 4i?i 3D) and = 
(/ix/?7io£sc)13.6eV, where mo is the free electron mass. 
The curves for the 2D and 3D binding energies are qual- 
itatively similar. However, the values in 2D are signif- 
icantly larger than those in 3D. This is physically ex- 
pected since the reduction of dimensionality allows for 
much stronger Coulomb interactions and more localized 
wave functions to enhance overlapping. We further no- 
tice that both, the 2D and 3D binding energies, have a 
minimum in the positronium limit, i.e., at mh/m e = 1; it 
then increases logarithmically as the mass ratio increases, 
until it saturates for large mass ratios. 

Our results gives a ground-state biexciton binding 
energy equal to 0.012i?^ D ' for mh/m e = 1 in 3D, 
which accounts for only about 40% of the more accu- 
rate variational result o 32 i 33 , while it reaches 0.21i?^? D ' ) 
when mh/m e — 1000, which accounts for about 70%. 
In 2D, our calculated binding energies give the ground 
state at 0.075i?^ 2D ' ) when m/,./m e = 1, which accounts 
for about 50% of the best variational result^, while it 
reaches 0.44i?^ 2D ' ) when ray l jm e = 1000, which accounts 
for about 80%. All this shows that our approach gives 
a much better ground-state energy when my l jm e ^> 1, 
i.e., close to the hydrogen molecule limit, possibly be- 
cause the exciton wave functions are less deformed when 
forming a molecule than in the case of lighter hole. 

One important advantage of the present approach over 
variational procedures is that it allows reaching the biex- 
citon bound and unbound excited states as easily as the 
ground state. Dashed curves in Figs. [7] and [5] show the 
binding energies of the biexciton bound states in 2D and 
3D. The number of bound states increases with the mass 
ratio mh/me, this number reducing to 1 for mh/m e < 20 
in 2D and mh/m e < 30 in 3D. For a large mass ratio 
mh/m e = 1000, we find 9 bound states in 2D and 8 
in 3D. The binding energy differences become smaller 
for higher excited states, evidencing a difference in the 
exciton-exciton interaction compared to the usual har- 




1 2 3 4 5 

p 



FIG. 9: (Color online) Plot of the ground state (rj — rjo) wave 
function L 2 |(r_i = 0, ri = 0, p|??)| 2 as a function of p in a^- 1 
unit, when the mass ratio mh/m e is equal to 5 and 50. Its p 
extension scales as the inverse of the biexciton Bohr radius, 
axx- We also show the same quantity for the first excited 
state (rj = r]i) when mh/m e — 50. 
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P 

FIG. 10: (Color online) Plot of ln|(r_i = 0, n = 0, p|ry)| 2 for 
three unbound biexciton states in a semilog plot, the mass 
ratio being nih/m e — 5 and the momentum p still in a^ 1 
unit. Note that the wave function is not here rescaled by L , 
so the amplitude of the wave function is significantly larger 
than for bound states, but its p extension significantly smaller 
for normalized functions. 



monic potential which leads to equal energy spacings be- 
tween cigenstates. 



B. Biexciton wave function 

In a previous work on biexciton, we have shown that 
the wave function of a biexciton made of opposite-spin 
carriers, with center-of-mass momentum K, relative mo- 
tion index 77, and electron and hole total spins S = 



10 



(S e ,S h ), splits as (sec Eq. (12) in Ref. M) 

(r ei ,r e2 ,r hl ,r h2 \K,r], S) = (R X x\K) (r_i, r x , u\rj, S), 

(40) 

where Rxx = (m e r ei + m e r e2 +m h r hl +m h r h . 2 )/2(m e + 
nth) is the biexciton center-of-mass coordinate, r x = 
r e 2 ~ r ?u an d r_i = r ei — r/j 2 are the electron-to-hole dis- 
tances of the two bright excitons having spins (±1), while 
u = (m e v e2 +m h v hl ) / (m e +m h )-(m e v ei +m h v h . 2 ) / (m e + 
rrih) is the distance between the center-of-masses of the 
two bright excitons. 

For bound biexciton, the wave function 
(r_i, ri, u|ry, S) has an extension of the order of 
the exciton size ax over r_i and over ri, and an 
extension of the order of the biexciton size ax x over u. 
By contrast, for unbound biexciton, the extension over u 
is as large as the sample size L, since unbound biexciton 
resembles very much an exciton with another free 
exciton moving around anywhere in the sample. Thus, 
in the case of bound states, dimensional arguments give 
the normalized relative motion wave functions through 

1 = J dridr_idu|(r_i,ri,u|?7, S)\ 2 



a x a x a xx \(0,0,0\v,SW 



(41) 



which leads to 



|(r_! =0,n =0,u = 0|t?,5)| . 

,a x axx 



D 



, (42) 



while, in the case of unbound states, axx is replaced by 



We can compute (i'i,i'j,p\rj,S) from (r_i, ri, p\r], S) 
through a double Fourier transform "in the exciton 
sense" (see also Eq. (31) of Ref. l3ll ). namely, 

(vi,Vj,p\v,S) = J dr_idri(i/i|r_i)(i/j|n)(r_i,ri,p|r/,5), 

(47) 

which also reads as 

(r_i,ri,p|r/,S) = ^ (r_i |^) (r x (v u v , p|r?, S). 

(48) 

Note that for r x or r_i equal to 0, the v exciton lev- 
els that survive in the above v sum are s-like states 
only. So, the (vi, Vj, p\r], S) function just is the function 

ftp (i>i,Vj) given in Eq. (|34[) . Since, when numerically 
solving the biexciton Schrodinger equation (|37l) for sin- 
glet state (S = 0), we have restricted the sum over v to 
the ground state v$, we must for consistency also keep vq 
only in the v sum of Eq. (|4"5)l . 

Figure H shows |(r_i = 0,r x = 0,p\v,S = 0)| 2 for the 
ground state when the mass ratio is mh/m e = 5, and for 
two bound states when the mass ratio is m/j/m e = 50. 
Note that Eq. (|45|) forces us to plot bound-state wave 



functions through L 2 |(r_i = 0,ri = 0, p|ry, S = 0)| 2 in 
order to have a quantity independent of sample size L. 
For unbound states, the |(r_i = 0,ri = 0, p|?y, S = 0)| 2 
function is peaked on momenta p^ which depend on the 
unbound biexciton energies. Figure [TU] shows that this 
peaked function broadens when the unbound biexciton 
energy increases, the broadening being due to exciton- 
exciton interactions. 



L; so, 



1^-1=0^1=0,^17,5)1- 



a\L 



(43) 



Biexciton absorption spectrum in pump-probe 
experiment 



To obtain (r_i,ri,p|?7, S), which is of physical rele- 
vance in photon absorption, we perform a Fourier trans- 
form as 

(r_i,n,p|77,S) = y<2u(p|u)(r_i,ri,u|r7,S), (44) 

where (p|u) = e' ip ' u /L D ^ 2 ; so, the extension over p of 
(r_i, ri, p|?7, S) is of the order ofl/axx for bound states 
(see Fig. [9]), and of the order ofl/L for unbound states 
(see Fig. nU)) . The same dimensional arguments then give, 
in the case of bound states, 



<r_! =0,n =0,p = 0| ? ;,S)| 2 - 



Wx L 



and, in the case of unbound states, 



|(r_! =0, n =0,p|r7, S')! 2 ~ ( -J- 

■ a x 



(45) 



(46) 



(i) Let us now first consider an initial state made of 
one circularly polarized photon cr + with momentum Q p h 
and frequency w, and one exciton already present in the 
sample, this (vq, Qi) exciton having an opposite circular 
polarization, cr_. After photon absorption, the final state 
contains two electron-hole pairs, their center-of-mass mo- 
mentum being K.; = Qpfc + Qj. The photocreated exciton 
interacts with the exciton present in the sample to pos- 
sibly form a biexciton. Since we are mainly interested in 
low-lying biexciton states, we shall focus on singlet states 
(S = 0). The Fermi golden rule gives the photon absorp- 
tion as (—2) times the imaginary part of the response 
function to one photon (w, Q„h) : This response function 
reads as (see Eq. (36) of Ref. 1311 ) 



Sxx (w, Qph', Qi) 

E 



uj + E, 



~g , (Qph+Qi) 



4AU 



(49) 



where Pi = (Qi — Q p / l )/2 is the relative motion momen- 
tum of the (X, X) pair and E Uo _Q i is the free exciton en- 
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crgy given by Eq. (|A.9[) . The biexciton oscillator strength 

fxx(p) in Ec i- <BS> is g ivcn b y 



/<&(p) = \il\ i L D \ y £(T = 0\ V )(u ,u,p\ V ) 



n\ 2 L D \(v ,r = 0,p|r?) 



(50) 



where f2 is the vacuum Rabi coupling, and Vj,p\rf) is 
the singlet biexciton wave function t/j, p\rj, 5 = 0). 
Using Eq. (l45l) . we can show that the biexciton oscil- 



lator strength /^jf (p) f° r P close to zero, is related 
to the exciton oscillator strength 

fx = M 2 L D \(r = 0|i,)| 2 ~ \n\ 2 (L/a x ) D (51) 
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FIG. 11: (Color online) Absorption amplitude A\^ X \uj,T) 



J XX 



0) 



defined in Eq. (|56[) . f° r the biexciton ground state in 2D, 



( a XX / L) D fx ■ (52) as a function of the photon energy u in R^ u ^ unit for var 

ious temperatures T. We find an asymmetric low-energy 
The prefactor (axx/L) D corresponds to the localization peak at -4.57i?^ D) , which is associated with the biexci- 



into a biexciton volume 



l xx 



the sample and initially delocalized over a volume L D . 
For very large sample size L, this a priori prevents using 
the same scale to draw bound and unbound biexciton 
absorption spectra through linear response to a photon 
field. 

(ii) We now turn to a more complicated situation in 
which one circularly polarized photon, er+, is absorbed 
in a dilute exciton gas having Nx excitons with opposite 
circular polarization, <r_ , which is what currently hap- 
pens in pump-probe experiments: One first prepares a di- 
lute exciton gas using a circularly polarized photon beam 
with low pump power; then we probe this gas with a 
weak photon beam having opposite polarization. We as- 
sume that these Nx excitons all are in the exciton ground 
state vq : the temperature being too small to have excited 
exciton states populated. For a dilute exciton gas, i.e., 
Nx (ax / L) D <C 1, we may neglect exciton many-body ef- 
fects, since these effects scale at least quadratically in the 
exciton density n x {= Nx/L D ). Indeed, the Pauli scat- 
tering for fcrmion exchanges between two excitons leads 
to terms in n 2 x . Such many-body effects would alter the 
absorption spectra presented below, because they affect 
the exciton energy states, and accordingly the biexciton 
state. In addition, the photocreated biexciton can in- 
teract with other excitons in the exciton gas. As a first 
approximation, we consider the dilute exciton g as as a 
set of noninteracting classical particles. The Qi exci- 
ton distribution for finite temperature T then is just the 
Boltzmann distribution 



of the exciton present in ton ground state having a binding energy 0.577^ D) below 



R 



(2D) 



4R^ D) . The hole-to-electron mass ratio is mj/m e 



5, as in usual GaAs samples, and 8tt\Q\ 2 N p hNx p/ L D — 1 is 
set equal to 1, the L D factor coming from the wave function 
part of this equation, due to Eq. (|45f) . 



a N(Qi, T) exciton distribution as 

S( xx Huj,T) = N ph J2N{Qi,T)Sxx(u,0;Qi), (54) 



with S X x(u, Q P h', Qi) given in ([49 



1. Bound biexciton 

By replacing the Q; sum in Eq. (|54[) by an integral 
with a constant density of states p for 2D systems, and 
by setting e = Q, 2 /2M X , we find that the absorption 
spectrum associated with a biexciton bound state ryo is 
given by 



4ir 2 \{l\ 2 N ph N x 
M x k B T 



pdee 



-e/k B T 



(v Q ,r = Q,p=y/M x e/2\r, ) 



xS(uj + s Uo -E va +e/2) 
For a photon detuning 5^ = £ Vo — e Va 



(55) 



leads to 



> 0, this 



N(Q h T) 



Nx_ 

L D 



2tt 



M x k B T 



D/2 



-Q 2 i /2M x k B T 



(53) 



normalized through J^q. N(Qi,T) = Nx- 

Since Q p h ~ on the characteristic electron scale, 
we can write the response function of N p h photons to 



,/yv,, „ 8TT 2 \Q\ 2 N v hNxP 9A<"0) II. T , . 

' Mxk b L 

(^r = 0,p^^MxS ( p l o) \ Vo ) 
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Figure [Til shows the absorption A\^ X ^ [uj, T) (|56|) asso- 
ciated with the biexciton ground state 770 in the presence 
of a dilute exciton gas as a function of the photon energy 
uj for various temperatures. We identify the biexciton 
ground state with the low-energy peak lying 0.57i?^ D '' — 
which is the biexciton binding energy — below a larger 
peak at — 4i?^ D ' ) , shown in Fig. [T2] that accounts for 
the photocreated exciton scattered by the free excitons 
present in the sample. Note the (L/a x ) D factor being 
the scale between the two figures. The lineshape of the 
ground-state biexciton peak is asymmetric with a long 
tail on the lower-energy side, the peak strength decreas- 
ing with increasing temperature. This low-energy tail 
reflects the fact that the photon energy required to meet 
the biexciton ground-state energy, as a result of energy 
conservation in Eq. (|55p , is smaller if the exciton already 
present in the sample has a larger kinetic energy. So, 
the biexciton binding energy corresponds to the upper 
sharp edge of the peak, since the tail comes from biexci- 
ton with larger relative motion momenta. Of course, this 
upper sharp edge is smoothed into a high-energy tail if 
we include a broadening of the delta function in Eq. (1551) , 
which amounts to replacing i0 + in Eq. ([4U]) by a finite 
width i'y. This low-energy tail is also found in the pho- 
toluminescence spectra of ID quantum wires 2 ^. 



2. Unbound biexciton 

Since the unbound biexciton energies are close to the 
exciton energies, photoabsorption with formation of an 
unbound biexciton A^ xx > is mixed with photoabsorption 
with formation of an exciton A^ x ^ , depending on the cap- 
ture rate / of an exciton by the photocreated exciton. As 
a first approximation, one can write the resulting absorp- 
tion spectrum as 

fA (xx ) + {1 _ f)A (x) = A (x )+f ^ A (xx) _ _ (57) 

For Nx excitons in a sample volume L D , the capture rate 
/ should be of the order of the exciton volume divided by 
the average volume occupied by one free exciton in the 
sample, namely, 

f^T^=N x m D . (58) 



L D /N X 

However, since excitons with different Qj's contribute 
differently to the biexciton absorption, as seen from 
Eq. ([54]). the absorption spectrum should in fact read, 
instead of Eq. (157 



/ ax 
\ L 



A(u,T) = A^ x \l 

x ~A( xx Xu,Qi 



Q, 

-A^(lu) 



(59) 



with N(Q it T) given by Eq. ([52]). We see that the ab- 
sorption spectrum reduces to the exciton spectrum in the 
absence of free excitons, Nx = 0, as physically required. 



X 

x 




3.97 



FIG. 12: (Color online) Absorption amplitude A (xx) (lu,T), 
denned in Eq. (|60l) . for the biexciton unbound states in 



2D, as a function of the photon energy to in i4 3D) unit for 
various temperatures. We find a large peak centered on 
We have taken an hole-to-electron mass 
ratio irih/m e = 5, a broadening 7 = 0.002i?^ D ' in Eq. (|61l) . 
and we have set 8n\Q\ 2 N p hNxp/a x = 1- 



-R 



(2D) 



-4R 



(3D) 



The exciton absorption spectrum A^ x >{oj) is made of 
delta peaks centered on the exciton energies eM, and 
weighted by the value at r = of the exciton wave func- 
tion squared |(r = 0|^)| 2 . When taking into account 
the finite exciton lifetime, these delta peaks broaden into 
Lorentzian functions with finite width. 

Although, according to Eq. (|59|) . photon absorption 
contains an exciton and a biexciton part, we have chosen 
to only show here the part of the spectrum coming from 
unbound biexcitons, namely, 



A^(u,T) = 



47T 2 \n\ 2 N ph N x 
M x k B T 



E j P det 



-s/k B T 



;i/ o ,r = 0,VAfe/2|»?) 



x8 (lu + e h 



f 



■12), 



(60) 



in order to avoid ambiguity coming from the relative 
weights of the exciton and biexciton lifetimes, which are 
sample dependent. To compute A^ xx \uj,T), we replace 
the discrete sum over r\ by a continuous sum over p,, and 
we introduce a finite lifetime by replacing the delta func- 
tion by a Lorentzian function having a small half-width 
7, namely, 



8(u>) 



uj 2 + 7 2 



(61) 



Figure [12] shows the lineshape of the peak located at 
— 4i?)^ D \ which is the Is exciton level in 2D quantum 
wells. This peak corresponds to unbound biexcitons, i.e., 
exciton-exciton scattering states. We see that this peak 
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spreads on both sides of the Is exciton level, due to en- 
ergy conservation enforced by the broadened delta func- 
tion. The peak lineshapc essentially is Lorentzian with a 
peak height slightly decreasing with temperature. This 
can be attributed to broadening at large relative motion 
momentum induced by exciton-exciton scatterings be- 
cause contributions from large relative momenta start to 
weight in when the temperature gets large. Due to energy 
conservation enforced by the delta function in Eq. (|60|) . 
the broadening induced by temperature may lead to un- 
bound states with a much broader lineshape on top of 
the sharp exciton peak. The broadening of the peak re- 
mains close to O.Oli?^ even though the temperature 

changes from O.Oliij,^ to 0.03i?^? D '. This indicates a 
"non-thermal broadening" behavior, since, for such tem- 
peratures, the intrinsic broadening due to exciton-exciton 
scattering is dominant over thermal broadening. 



VI. CONCLUSIONS 

We have constructed the Schrodingcr equation for 
biexciton written in the exciton basis, instead of the stan- 
dard free-carrier basis. To choose such a basis is moti- 
vated by the physical fact that Coulomb interaction be- 
tween electron and hole is far stronger than between two 
dipole-like excitons. However, in this exciton formula- 
tion, we must handle the fact that two excitons can ex- 
change their carriers. The exact handling of this carrier 
exchange, as induced by the Pauli exclusion principle, is 
made possible thanks to the recently developed compos- 
ite boson many-body formalism. 

By restricting the exciton levels to the ground state, 
we have numerically solved the resulting biexciton 
Schrodingcr equation for the ground state, as well as for 
the bound and unbound excited states in 2D and 3D sys- 
tems. Our results for the ground-state binding energies 
agree reasonably well with those obtained from varia- 
tional methods. The main advantage of the present ap- 
proach is to easily reach excited states, which are out of 
reach from usual variational procedures. We then use the 
obtained biexciton wave functions to calculate the optical 
absorption spectrum in the presence of a dilute exciton 
gas in 2D. The spectrum shows a small asymmetric low- 
energy peak associated with the biexciton ground state, 
and larger peaks coming from exciton-exciton scattering 
states. 
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Appendix I. EXCITON OPERATOR IN 
MOMENTUM SPACE 

The index i in a state \i) of an exciton made of electron- 
hole pairs in a translationally invariant system, refers to 
the pair center-of-mass momentum Qi and the internal 
motion index of the exciton state at hand, Ui, i.e., i = 
(vi> Qi)- The wave function of this (i>i, Qj) exciton splits 
as 

(f>i{r e ,r h ) = (R|Qi) (r|^), (A.l) 

where R = (m e r e + m^r^)/ (m e + m^) is the center- 
of-mass coordinate of the pair and r = r e — th is the 
distance between electron and hole. The center-of-mass 
momentum wave function is a plane wave, (R|Qj) = 
e iQ f R^D/2 f or a g j zc £ sample. By writing the rela- 
tive motion wave function as (r\v) = ■ ( r \Pi) (PiW) ■, 
we can rewrite <fii(r e , r/j) as 

Mre,r h ) = 2_>ik> — 7^72 7-572 > 

(A.2) 

where we have set 

m e 

a e = = 1 - a h . (A. 3) 

m e + mu 

So, the electron and hole momenta read in terms of the 
center-of-mass and relative motion momenta of the pair, 
Qi and p^, as 

k e = Pi + a e Qi, k h = -pi + ahQi- (A. 4) 

The exciton relative motion wave function (k|i/j) obeys 
the Schrodingcr equation 




where [i x x = m e 1 +m h 1 is the electron- hole pair reduced 
mass, while the Coulomb potential V^j is given by 

J 2ire 2 /e sc L 2 q in 2D 
V «~ \ 47re 2 /e sc LV in 3D ' [A - b > 

From Eq. (|A~2|) . we can deduce the creation operator 
of the i exciton as 

B t = 4 ilQi = E^i^U.qA+^q,; ( a - 7 ) 

k 

This operator creates a single-pair eigenstate 

(H-Ei)Bl\v) = 0. (A.8) 

The exciton energy Ei contains a relative motion part 
and a center-of-mass kinetic energy, namely, 

Ei=E Vi , Qi =s Ui+ ^-, (A.9) 
where Mx = m e + rrih- 
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Appendix II. ON THE VARIOUS 
SCATTERINGS OF TWO EXCITONS 



In this appendix, we rederive various scatterings ap- 
pearing in the biexciton Schrodinger equatio n 22,29 . With- 
out any loss of generality in calculating these scatterings, 
we set the biexciton center-of-mass momentum K to be 
zero. The two excitons then have opposite momenta. 
This just amounts to working in the reference frame of 
the center-of-mass of the exciton pair. 



A. Pauli scattering 

Let us first consider the Pauli scattering for hole ex- 
change between excitons starting in "in" states and 
ending in "out" states (m, n), the excitons m and i having 



the second diagram of Fig. |3l it reads as 



the same electron, as shown in the diagram of Fig. 1(a) 
This scattering is given, as read from the figure, by 

A » (Zi) =EHkX)(n|p^}(p/ l Pe|i)(k ft k e |i} 1 

(B.l) 

Similarly, the Pauli scattering for electron exchange 
shown in the diagram of Fig. 1(b) | with the excitons m 
and i having the same hole, is given by 

X * (ml) = EHPX>MW)<P*Peb'><kAkel*'>- 

(B.2) 

It is easy to check that A e ( ™ \ ) = \h(™l), the Pauli 
scattering A^ ( n 3 . ) being often written as A ( n 3 . ) for 
simplicity. 

For two excitons i and j having opposite momenta, i.e., 
i = Q) and j = {vj, — Q), Eq. (jB.ll) reduces to 



^^ m |k')(^| P ')(p|^)(k|- 



(B.3) 



provided that the (k',p',k, p) momenta are such that 
-p' - a h Q! = -k + a h Q, -k' + a h Q' = p a h Q, 
k' + a e Q' = k + a e Q, and p' — a e Q' = p a e Q, as read 
from Fig. [2] When inserted into Eq. (jB.ll) , we get 



A; 



( 0"»,-Q') in, 



-Q) 
,Q) 



P P 

x(P+-^-|^)(p- -g^n 



where the momenta P± are defined by P± 
Q')±a e (Q'-Q). 



(B.4) 



B. Direct Coulomb scattering 

We now turn to the part of the direct Coulomb scat- 
tering resulting from hole-hole interaction. As seen from 



i-dir I n j \ _ 
^hxh 2 \ m i ) 



E V-qHke.kft-q) (B.5) 

q,ke,kfe,p e ,Ph 

x(n\p e ,p h + q)(p ft ,p e |j)(k h ,k e |i). 



In this direct process, the excitons keep their electron and 
hole components. As seen from Fig. [U for two excitons 
i — (vi, Q) and j — (uj, — Q), Eq. (|B.5|) reduces to 

^y q ^ m |k')K| P ')(p|^)(k|^), (B.6) 

provided that the (k',p',k, p) momenta are such that 
-k'+auQ' = -k+a^Q-q, -p'-auQ' = p a h Q+q, 
k' + a e Q' = k + a e Q, and p' — a e Q' = p a e Q, as read 
from the figure. When inserted into Eq. (|B.5|) . this gives 



tdir / (f„,-Q') (t^,-Q) \ CR 71 

<*hih 2 \ (i/ m ,Q') (yi,Q) J v '> 

= Vp E( I/ m|k - a e Po)(^n|p + a e Po)(k|fi)(p|i/j), 



k P 



where Po = Q' — Q is the momentum transfer. Using the 
same procedure for £eje 2 > £fc ir e , an d ^t^k 2 ' we enc ^ w ^ n 
a direct Coulomb scattering between two excitons given 
by 

tdir I (fn,-Q') (^j,-Q)\ (td Q\ 

« V C"m,Q') ("<.Q) J [ ' 

= Vp E [(^m|k + a, l P )(^, i |p - a h P ) 
k P 

+ (i/ m |k - a e P )(v n \p + a e Po) 
- (v m \k + a h P ) (v n \p + a e P ) 

-{v m \k - a e P a )(is n \p - a h P ) 

x(k|^)(p|^-). 

By noting that 

E[^'|k + a/»q>-(^|k-« e q>]<k|f> 

k 

= {y'\e ia ^ T - e- ia «* T \v) = T v -M), (B.9) 
this direct Coulomb scattering splits as 

^ dir ( S7q?) = ^-^(Po)r^(-Po). 

(B.10) 

For and i/ restricted to the exciton ground state vq, we 
find 



/a h a x q\ ( 

> (q)= n^— )- g { 



OeOxq_\ 

2 J ' 



(B.ll) 



where, for 2D and 3D systems, g2D(p) = (1 + p 2 /4)~ 3 / 2 
and <73d(p) = (1 +p 2 )~ 2 - We see that 7I/ oi , (q) depends 
on the magnitude of the momentum transfer q = |q| only 
and that 7^ 01 , (q = 0) = 0. We also note that % QV0 (q) = 
for a e = ah, i.e., equal electron and hole masses. 
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C. Exchange-Coulomb scatterings 

Two excitons can also have exchange-Coulomb scatter- 
ings. These are defined as 



r(:A) = EM^K dir Ci), (B.i2) 

rs 
rs 

depending on if the exchange takes place after or before 
Coulomb interaction. 

The part of the "in" exchange- Coulomb scattering due 
to hole-hole interaction, as shown in the first diagram of 
Fig. [5J is given by 

tthAZi) = E V n (m\k e ,p h + q) (B.14) 

q,k c ,k h ,p e ,pfc 

x(n\p e ,k h - q)(p /l ,p e |j)(k /l ,k e |i). 

As seen from Fig. El for i = (z/j, Q) and j = (i>j, — Q), 
Eq. (|B.14[) reduces to the same equation as Eq. (|B.6|) , the 
(k',p',k, p) momenta being now such that k' + a^Q' = 
p - cthQ + q, p' - a h Q' = k + a h Q, - q, -k' + a e Q' = 
p — a e Qj an( l p' — a e Q' = k + a e Q. This leads to 



An ( (f„,-Q) {uj, -Q 



Q) (f«,Q') 

£ y p (^jk+^±P)K|k-^±P; 

k,p#0 



x(k+ |^)(k |i/ 



(B.15) 



with the lower sign in front of p for c\C2 = /lift. 2 and the 
upper sign for C1C2 = e.\&2- 



The same procedure for the electron-hole part of the 
interaction yields 



tin / (f»,-Q') ("j, 
y (v m ,Q,') (vi 



-Q) 
,Q) 



- 2 F p (^ m |k+ P ^^)K|k- P ^-tP) 



k,p#0 

x(k+ |^}(k 



(B.16) 



with the upper sign for c\d,2 = ei/12 and the lower sign 
for C\d2 — h\e2- 

Note that we can eliminate the sum over p in the 
electron-hole part of the £j, n d scattering, by setting k' = 
k p/2 and by using Eq. We then find 



An ( {u n , 

= E 



-Q ) (fi ,-Q) 
,Q') (^,Q) 

(k + P_/2) 2 \ ; „ P_. 



xK|k-^)(k+^t|^)(k-^t|^-){B.17) 
Similarly, by setting k' = k + p/2, we find 



tin / {v n 



-Q') (^,-Q) 

l,Q') ("i,Q) 



(k-P_/2) 2 \ 7 P_ 

2^ I £ ^ - I ^™l k + — 



2fix 



■> 



x<i/„|k- ^)(k+ ^h)(k- ^h)(B.18) 
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